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ABSTRACT 

We present a new sample of 4634 eclipsing binary stars in the Large Magellanic 
Cloud (LMC), expanding on a previous sample of 611 objects and a new sample of 1509 
eclipsing binary stars in the Small Magellanic Cloud (SMC), that were identified in the 
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light curve database of the MACHO project. We perform a cross correlation with the 
OGLE-II LMC sample, finding 1236 matches. A cross correlation with the OGLE-II 
SMC sample finds 698 matches. We then compare the LMC subsamples corresponding 
to center and the periphery of the LMC and find only minor differences between the two 
populations. These samples are sufficiently large and complete that statistical studies 
of the binary star populations are possible. 

Subject headings: binaries: eclipsing — Magellanic Clouds — surveys 



Introduction 



Eclipsing binary stars (EEs) are important for astrophysical research in man y ways. They 



may be used to obtain accurate estimates of star masses and radii (I AndersenI 1 1 99 ll . and references 
therein). Precise determination of stellar parameters can in turn be used to put theories of stellar 
structure a nd evolution to a stringent test by comparing measure d parameters with theoretical 



predictions (jLastennet Sz Valls-Gabaud 



2002 



Lastennet et al.ll2003l . and references therein). 



EEs may also be use d for distance determina t ion a nd this use goes back se veral decades; 
its history is reviewed by iKruszewski &: SemeniukI (jl999l ) . Since IStebbind (Il91lh used an esti- 
mate of the parallax to /3 Aurigae to infer the surface brightness of both its components, it has 
been known that a good photometric light curve plus a double line spectroscopic orbit admits a 
simple g eometric r elatio nship between the surface brightnesses of the stars and the distance to 
the EE; IStebbind (jl91ll ) however had no way at the time to make the reverse "surface bright- 



ness to distance" i nference and his pa.per does no mention this po s sibility. Afte r 



other early papers (jCaposchkin 



1933 



Woollev 



1934 



Pilowski 



1936 



Kopallll939l : iGaposchkin 



Stebbini 



(1191 ih 



193S 



19401') used parallaxes obtained in dependently to estimate surface brightnesses, but, as remarked 
by IKruszewski &: SemeniukI (jl999l ) , these pioneers surely knew of the potential of this te chnique to 
estima te distances. Modern analyses of EEs have usually focussed on this technique (e.g. lAndersen 
(jl99ll )). The method affords high precision due to its purely geometrical nature and has been ap- 
plied by a number of authors t o determine the distance to the Large Magellanic Cloud (LMC) usiii g 
H V2274 dudalski et al.lll998al: Guinan et al.lll998l:lNelson et aklboool : ICroenewegen fc SalarijboOllv 
HV9 82 (|Fizpatrick et ahlbood ). EROS 1044 teibas et aklbood ) and HV 5936 dFitzpatrick et. al 



2003); an attempt to use EEs to determine the distance to M31 is currently under way ( Ribas et al 



2003t ) and the DIRECT project is attempting to ineasur e the distance t o M31 and M33 via EEs 
and Ce pheids (iKaluznv et al.lll998 : lEonanos et al.l l2003l : lEonanosI l2005l ) : other recent examples 
include iMichalska &: Pigulskil (|2005|) who pre sent a sample of detached binaries in the LMC for 
distance determination and lRibas et al. (2005) who present the first determination of the distance 
and properties of an EE in M31; iNorthI (|2006l ) presents a sample of EEs with total eclipses in the 
LMC suitable for spectroscopic studies. In general it is important that distances be determined 
using a large sample of EEs to minimize the impact of systematic errors. A recent collection of 
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references on extragalactic binaries can be found in lRibas &: Gimenea (|2004l ) 



Large-scale surveys to detect gravitational microlensing events have identified and collected 
light curves for large numbers of variable stars in the bulge of the Milky Way and in the Magellanic 
Clouds. Eclipsing binary stars comprise a significant fraction of these collections. The MACHO 
collabo ratioru has presente d a sample of 611 EEs in the LMC with preliminary analyses of their 
orbits (jAlcock et a l]|l997a|) . A catalogue of 3031 EEs in the LMC found in the MACHO database 
has been just published bv lDerekas. Kiss. &: Eeddind (j2007l ): this catalogue was compiled by ana- 
lyzing a list of 6835 stars classified as possible EEs in the MACHO database; a cross correlation 
between these 6835 stars and our sample finds just 1987 matches, thus at least about 2700 EEs in 
our catalogue are new identifications. The 6835 classified as possible EEs were found in regions of 
parameter space such as color, magnitude, and period, where one does not expect to find pulsating 
variables and therefore the detected variability of these stars was tentatively ascribed to eclipses. 
Regions where pulsating variables could exist were not considered while making this preliminary 
classification and EEs there were therefore not included in the list. In our search we did not rely 
primarily on cuts in parameter space and and we did not exclude a priori regions of this space 
where pulsating variables are present; therefore we were able to classify many EEs in these regions, 
that were not included in the^reU mmarv dasdflcation. The OGLE c„llab„™lio,a has introduced 
a sample of 2580 EE s in the LMC dWvrzvkowski et al.l l2003l') and of 1351 EEs in the Smah Mag- 
ellanic Cloud (SMC: IWvrzvkowski et al.l 120041') . Eoth samples were selected from their catalogue 
of variable stars in the Magellanic Clouds (IZebruh et al.l 120011) compiled from obser vations taken 



via Difference Image Analysis (DIA: IZebruh. Soszyhski. &: Wozniak 



during the second part of the projec t (OGLE ILlUdalski. Kubiak &: Szvmahskilll997l l and reduced 



200ll)^ An earlier samp le of 79 

EEs in the bar of the LMC was presented by the EROS collaboratioiPl Grison et al. 1995 ). Other 
large variable star data sets are being produced by surveys not specific ally designed to detect 
gravitational microlensing, such as the All Sky Automated Survey (AS AS: IPoimahskil[l997l ) 0. 



The availability of large samples of EEs (and the even larger ones that can be found by future 
surveys such as Pan-STARRd3 and LSSTO) can have an important impact on stellar astrophysics. 
This impact can arise in two qualitatively different approaches. First, a large catalogue allows 
the discerning researcher to select carefully a few EEs for detailed follow-up study; the distance 
estimation described above is an example of this. Second, statistical analyses of an entire population 
become possible when a large collection is assembled; such analyses of EEs have not previously been 



77/www . macho . mcmaster . ca/ 1 
^http: //sirius . astrouw. edu.pl/~{}ogle/| 
Ihttp : //eros . in2p3 . f r/ 1 

' http: //www. astrouw. edu.pl/~{}gp/asas/asas .html 
77/paii- Starrs . if a.hawaii . edu/public/| 



* ]http: //www. Isst . org/lsst_home . shtml/| 
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possible. To fulfill this promise there are challenges to overcome, including finding EBs in large 
data sets and automating their analysis. With regard to the first task, the discovery problem is 
complicated by the fact that EBs do not have clear relationships between their parameters (period, 
luminosity, colors) as do the major classes of pulsating variables. This makes it difficult to find 
them via simple and well understood cuts in parameter space. The first step toward automated 
discovery is thus to have a large sample of data on which to experiment with search techniques. 
This non-trivial exercise in mining large data sets can be useful for future surveys t h at ar e not 
necessarily aimed a t bina ry star research. An example is given by IWyrzykowski et al.l (120031 ) and 



Wvrzvkowski et al.l (|2004l ) who employ an artificial neural network to identify EBs in the OGLE-II 



LMC and SMC samples, but more needs to be done. With regard to analysis of EBs, the traditional 
approach has been to carefully analyze indi vidual systems with the he l p of dedicated computer codes 



such as the Wilson-Devinney code (WD: IWilson Devinnev 



1971 



Wilsonlll979l ). This becomes 



impracticable when many thousands of stars are involved and an automated approa ch is required. 



The li ght curves in a previous sample of 1459 EBs in the SMC foun d by OGLE-II ([Udalski et al 



1998bl ) were systematically solved by IWyithe &: WilsonI (|200ll . l2002l ) using an automated version 



of the WD code; the ASAS collaboration has developed an aut o mated classif ication algorithm for 
variable stars based on Fourier decomposition (|Pojmahskill2002l ): iDevoii (120051 ) found and analyzed 
10000 Bulge EBs from OGLE-II using DEBilfl, an EB analysis code that allows automated solutions 
of large EB data sets and works best for detache d EBs; a genetic algorithm based approach to finding 
good initial parameters for WD is described in iMetcalfd (|l999l ). 



This paper is the first of a series of papers aimed at describing the EB samples in the MACHO 
database and is organized as follows: Section [2] introduces the LMC and SMC samples; Section 
[3] describes the Color Magnitude Diagram (CMD) and the Color Period Diagram, pointing out 
significant features in them, Section [J] compares the LMC and SMC samples, Section [5] describes 
the results of the cross correlation with the OGLE LMC and SMC samples, and Section [6] reports 
where and in what form the data presented in the paper can be accessed on line. 



2. The Samples 

2.1. The MACHO Project 

The MACHO Project was an astronomical survey whose primary aim was to detect gravita- 
tional microlensing events of background sources by compact objects in the halo of the Milky Way. 
The gravitational background sources were located in the LMC, SMC a nd the bulge of the M ilky 



Way; more details on the detection of microlensing events can be found in lAlcock et al.l (l2000al ) and 
references therein. Observations were carried out from July 1992 to December 1999 with the dedi- 
cated 1.27m telescope of Mount Stromlo, Australia, using a 2 x 2 mosaic of 2048 x 2048 CCD in two 



http: //www. cf a.harvcird. edu/~{} jdevor/DEBiL.html 
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bandpasses simultaneously. These are called MACHO "blue", hereafter indicated with 144achO) 
with a bandpass of ~ 440 — 590nm and MACHO "red" , hereafter indicated with -Rmacho j with a 
bandpass of ~ 590 — 780nm; th ese widths are between the half-response points as estimated from 
Figure 1 of (jAlcock et al.lll999l ). The bandpasse s and the transform ations to standard Johnson V 
and Cousins R bands are described in detail in (lAlcock et al.lll999l ): see in particular their Figure 
1 for the instrumental throughput of the two MACHO bands. Each MACHO object is identified 
by its field number (1 — 82 for the LMC, 201 — 213 for the SMC), its tile number (which can over- 
lap more than one field), and its sequence number in the tile. These form the so called MACHO 
Field. Tile. Sequence (FTS), which is used in this paper to label EEs. Note that, since some overlap 
exists between fields, one star may have two or more FTS identifiers. 



2.2. The Large Magellanic Cloud Sample 

The LMC sample we present comprises 4634 EEs selected by a v ariety of methods wh ich we 



describe in this section; the sample includes the 611 EEs described in lAlcock et al.l (|l997al ). The 



LMC magnitudes quoted in this paper have been obtained by using the following transformation: 
V 



R 



HdACHO + 24.22 mag - 0.18(yMACHO - -Rmacho) 
-Rmacho + 23.98 mag + 0.18(Vmacho - ^macho)- (1) 



From now on we will use the symbols V, R, and V — R to refer to standard magnitudes obtained 
from Fmacho and -Rmacho via Eq. [1] for the LMC and Eq. [2] for the SMC, and not corrected 
for reddening; we will also use Vmacho and -Rmacho to indicate instrumental magnitudes. The 
observations number in several hundreds in both bandpasses for most light curves; Figure [1] shows 
histograms of the number of light curve points of the EEs in both bands; the I^^ACHO band has on 
average more observations than the -Rmacho band because one half of one of the red CCDs was 
out of commission during part of the project. The central fields of the LMC were observed more 
often and the periphery less often as shown by the three peaks in the distribution where the first 
peak corresponds to the LMC periphery and the other two correspond to the center. 



2.3. Identifying Eclipsing Binary Stars in the MACHO database 

This Section describes the techniques employed to identify EEs both in the LMC and in the 
SMC; the results we quote are relative to the LMC. From now on we will always use the term 
unfolded light curve to indicate a set of time ordered observations and will reserve the term light 
curve to indicate a set of time ordered observations folded (or phased) around a period, omitting 
for brevity the adjectives "folded" and "phased"; we will also use the terms EE, system, and object 
interchangeably. 



All sources in the survey were subjected to a test for variability (jCook et al 



19951 ) and a large 
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Fig. 1. — Histograms of the number of light curve points for both bands for the LMC sample. The 
first peak corresponds to objects in the LMC periphery which was observed less often than the 
center, the other two peaks to the regions in and near the central bar. 



-7- 



number of variable sources were identified (|Alcock at alJll995l : lAkock et alJEoOGal lbl. Il997bl l. This 



first test starts by first eliminating the 20% most extreme photometric data points; the resulting 
unfolded light curve is fitted to constant brightness and its x^/dof is calculated. The elimination 
of the most extreme points is expected to reduce the influence of noise and yield a good fit for a 
constant source, but not for a truly variable one. We consider the source variable if the x^/dof thus 
computed can occur by chance with probability 1% or less. Sources that were flagge d for variability 



were tested for per iodicity. P e riods were found using the Supersmoother algorithm (lReimannlll99^ 



first published by iFriedmanI (jl984l )). The algorithm folds the unfolded light curve around trial 
periods and selects those periods in which the smoothed light curve matches the data best in a 
statistical sense; we have selected the best 15 possible periods ranked by the smoothness of the 
light curve. Periods were found separately for the red and blue unfolded light curves. The period 
selected as the best one by the program turned out to be "correct" in 88% of the EBs for at least 
one band; for 10% of the EBs the second best period turned out to be correct for at least one 
band and only for < 0.5% of the EBs did the procedure fail to find a good period. In these cases 
"correctness" was determined by direct visual inspection. 

The Supersmoother program can fail in two manners when fitting EBs. First when one eclipse 
(the secondary) is very shallow, Supersmoother may not recognize it and yield a period twice the 
correct one. Second, when the two eclipses have nearly equal depth Supersmoother may confuse 
the secondary and the primary eclipses yielding a period half the correct one. The first failure 
happened, for one or both bands, in about 2% of EBs, whereas the second happened in about 13% 
of EBs. These cases are easy to correct upon visual inspection. For ~ 100 EBs Supersmoother gave 
a period which was some other multiple of the correct one for at least one band; these were fixed 
upon visual inspection. For the remaining EBs we tried folding the light curves around the other 
periods selected by Supersmoother and managed to identify the correct period for most of them. 
In 27 cases in which there were OGLE-II counterparts we adopted OGLE periods since, though 
differing in some cases by less than 1% from the periods found by Supersmoother, they gave a 
much better light curve. We found 51 stars in which the secondary eclipse was not evident, either 
because it was shallow or because the light curve was noisy, but with an OGLE-II counterpart in 

Table 1. Period determination via Supersmoother. 



Color 


Number of EBs 


= Pss* 


= 2Pss* 


Pt = Pss/2t 


Other 


Red 


4634 


3432 


604 


103 


495 






74% 


13% 


2% 


11% 


Blue 


4634 


3460 


583 


108 


483 






75% 


13% 


2% 


10% 



^Real period. 

*Best period found by Supersmoother. 
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which it was clearly visible; these stars have not been included in the catalogue. The periods in 
the two bands differ on average by 0.02%. These results are summarized in Table [TJ 

The search for variable objects in the LMC gave ~ 207, 000 objects of which ~ 66, 000 were 
found to be periodic. To find EBs in this sample we considered a variety of properties of light 
curves. The techniques we employed are: 

1. Look at the number of photometric excursions ("dips") in the light curve. An EB is expected 
to show two "dips" in an entire period corresponding to the two eclipses as opposed to a 
Cepheid or an RRLyrae star for which only one is expected. The number of photometric 
excursions was calculated by Supersmoother by counting the number of times the smoothed 
light curve crosses the mean: we selected stars with two excursions in both bands. We addi- 
tionally imposed a cut on light curve amplitudes: calling Amply and AmplR the amplitudes 
of the blue and red light curves respectively, as computed by Supersmoother, we imposed 
Amply /AmplR < 1.2. This amplitude cut was imposed to help in eliminating RRLyr^s from 
the sample, since, considering a population of several thousands probable RRLyr^s found in 
the MACHO database, we found that, on average, Amply /Amplj^ 1.27; imposing a cut 
Amply/AmplR < 1.2 should therefore filter out many RRLyrae. In our sample just 2421 
EBs (~ 52%) pass this cut; we then removed the amplitude cut and look at the number of 
photometric excursions alone we found that 3039 EBs (~ 66%) pass this relaxed cut. 

2. Look at the ratio of number of points five standard deviations away (s5) from the median to 
the number of points five standard deviations below the median (s5d) . This ratio is expected 
to be ~ 2 for a typical single variable star. For an EB we expect this ratio to approach ~ 1 as 
the signal to noise in the photometry increases, as most "outlier" points are due to eclipses. 
We imposed a cut s5/s5d < 1.2 in both bands and found that 3417 EBs pass it (~ 74%), 
whereas for the overall variable star dataset the figure is ~ 19, 600 out of 66, 000 or 30%. 



3. Use a decision tree. We applied the decision tree program described in lMurthv. Kasif &: Salzberg 



(| 19941 ). which was run on all the variable objects on the catalogue and gave for each the prob- 
ability that it was an EB, an RRLyr^, a Cepheid, a long period variable or an unknown 
object. In all ~ 17000 objects were found most likely by the decision tree to be EBs but only 
3281 were found to be real on visual inspection and included in the sample. 

4. Use a similarity t echni que. We finally tested the sample with a technique described in 



Protopapas et al.l (|2006l ) aimed at finding "outliers" in large data sets of variable star light 
curves. The technique aims at finding objects whose light curves are most dissimilar, in a 
statistical sense, from an "average" light curve built out of all the light curve in the data set. 
This is accoi nplished by lookiiig at a. 11 the pairs of light curves to find their mutual similarity 



as defined in lProtopapas et al.l (120061 ) and, for each light curve, by then combining these mea- 
sures, to find its overall similarity to the rest of the sample. Light curves with low measure 
of similarity are fiagged as outliers. This approach was useful in finding misfolded lightcurves 
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since it found many objects for which the period for one band gave a badly folded light curve 
but the period for the other band gave a good folding. This happened for 198 EBs, despite 
these periods differing on average by just 0.18%; in this case we selected the period that gave 
the good light curve for both bands. 

Since all the techniques we used give some false positives, each candidate was also visually 
inspected before inclusion in the sample; we paid closer attention to those stars which could more 
easily be classified as EBs without being so, like ellipsoidal variables (see Subsection 12. 6p and 
Cepheids and RRLyrae mistakenly folded around a period twice their real one. 

These results of our search techniques are summarized in Table [21 and Table [3] shows the 
number of EBs that pass more than one cut. As the numbers show these heuristic tests are far 
from perfect and tend to give too many candidates; however we feel that the large dimension of our 
sample can allow the determination of more stringent tests, an absolute necessity for the analysis 
of future surveys. 

Our search gave 266 objects observed in more than one field: these duplicates have different 
MACHO field numbers, but typically the same tile number. In this case we summed the numbers 
of observations in both bands for each field and chose the one which had the highest total number 
of observations: the object is identified by that corresponding FTS only. 

Figure [2] shows a logarithmic histogram of the period distribution. Note that the periods range 
from a fraction of a day to several hundreds of days. Figure [3] shows a histogram of the distribution 
of median V, R and V — R. Note that the magnitudes range in values from ~ 20 mag to ~ 14 mag 
both in V and R with a peak around 18 mag. The average photometric error for the LMC is 
~ 0.05 mag in both instrumental bands; the average error for a light curve as a function of median 
relative magnitude is in standard magnitudes is shown in Figure [H 

2.4. Root Mean Square of residuals for the LMC sample 

For the LMC sample we estimated the distribution the Root Mean Square (RMS) of the 
residuals Oi — Ci of the observations around a theoretical light curve as a function of median 



Table 2. Summary of cuts applied to define the LMC sample. 





Number of EBs 


2 dips 


2 dips and Amplv < 1.2AmplR 


s5/s5d < 1.2 


Decision tree 


EB 


4634 


3039 


2421 


3417 


3281 






66% 


52% 


74% 


71% 


All variable sources 


~ 66, 000 


~ 8000 


~ 7300 


- 19, 600 


- 7900 






~ 12% 


~ 11% 


-30% 


~ 12% 
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Fig. 2. — Period histogram for 4634 EBs in the LMC sample; the distribution peaks strongly in 
the 0.8 — 4d range and has a tail in the 10 — lOOd range. The size of the bins is ~ 1/100 of the 
span of the logarithms of the periods. 
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Fig. 3. — Standard V and R magnitudes and V—R histograms for the LMC sample; the bin size 
is 0.1 mag for the V and R histograms and 0.01 mag for the V — R one. The V — R histogram is 
strongly peaked around V — R ^ mag, showing a majority of unevolved EBs, but the longer tail, 
with the smaller bump around V — R ^ 0.5 mag shows a sizeable minority of fairly evolved systems 
as also shown by the higher V histogram values around 16 mag with respect to the R histogram. 
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Fig. 4.— 



Errors as a function of median magnitude for the LMC sample. 
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relative magnitude where Oi is the value of the observed magnitude orbital phase (pi , and is the 
theoretical value at the sai ne phase. Theoretical values were obtained by fitting the light curves 



using the JKTEBOF@ code (|Southworth. Maxted fc Smal levll2004al: ISouthworth et al.l[2004bl l. The 



JKTEBOP co de is based on th e EBO P code (|Etzel.l981 : iPopoer k Etzel 



1981 



which implements 

the model by iNelson &: DavisI (|l972l ) with some modifications; JKTEBOP in turn adds several 
modifications and extensions to the original EBOP code that make it easier to use, especially when 
fitting a large number of light curves. Before fitting we eliminated outlying points by taking averages 
of all points in boxes containing from ~ 10 to ~ 20 points along a light curve and discarding the 
points more than 2 standard deviations away from these averag es. We obtained starting values for 
the model parameters by running the DEBiL code (|Devoiil2005l ) a nd using t he values it computed; 
the limb darkening values for the V and R bands were taken from (jCoxll2000l ) . We fixed the ratio of 
the masses, q, to 1, 0.1 and 10 and did the fit in each case taking in the end the best result. Finally 
we selected light curves with x^/dof < 2 to show in Figure [5j Out of 4636 EBs in the LMC the 
program converged in 4090 cases in the R band and in 4312 cases in the V band; we found 3067 fits 
with x^/dof < 2 in the R band and 3198 in the V band. We point out that those fits were made 
only with the aim of obtaining a good theoretical light curve for as many observed light curves as 
possible in a fast and automated manner, so that a residuals distribution could be calculated. In 
particular we did not attempt to accurately determine astrophysical parameters for our EBs. This 
is also the reason why we discarded points at just 2 standard deviations away from the moving 
averages, which could result in eliminating potentially interesting information for some EBs; such 
objects are obviously deserving of more in depth study which we did not attempt here. While in 
general our fits were good in the case of largely separated, undistorted systems, they were often 
bad for close, strongly distorted ones, which is to be expected since JKTEBOP is not meant to be 
used for such systems; also in several cases our fits were bad because the scatter of the observed 
values was larger than the observational errors, which suggests that other physical phenomena, 
such as pulsation of one or both components, are present. The RMS distributions of the residuals 



' jhttp: //www, astro .keele . ac .uk/~{}jkt /codes/ jktebop. html] 



Table 3. Summary of EBs passing more than one cut. 



Cuts 



Number of EBs 



2 dips and s5/s5d < 1.2 2061 

2 dips, Amplv < 1.2AmplR, and s5/s5d < 1.2 1686 

2 dips and decision tree 2031 

s5/s5d < 1.2, and decision tree 2574 

2 dips, Amplv < 1.2AmplR,, and decision tree 1761 

2 dips, s5/s5d < 1.2, and decision tree 1488 

2 dips, Amplv < 1.2AmplR,, s5/s5d < 1.2, and decision tree 1280 
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vs. median magnitudes V and R are shown in Figure [5j 

2.5. Examples of light curves 

Figures [U [71 EJ [U [lOl and [H] show some examples of hght curves. The panels on the left show 
the original light curves, those on the right show the light curves with the outlying points removed, 
the error bars, the theoretical light curves from the fit and the residuals The EBs shown are meant 
to be representative of the sample, this is why some examples of bad fits are included. Comparing 
the panels on the left with the panels on the right for these figures gives an idea of the effect of 
removing outlying points. In particular the figures suggest that for those objects with good fits 
the procedure resulted in the elimination of truly outlying points; these EBs are mostly detached 
with undistorted components. For objects with bad fits, which mostly comprise EBs with close 
and strongly distorted components, the situation is less clear. For example the system labelled 
1.3442.172 shown in Figure [7] exhibits some points in the R band, at secondary eclipse around 
phase ~ 0.4, that run almost parallel to the main light curve but at a higher magnitude. Such 
points may or may not be physically significant; some of these are removed by our procedure, but 
the fit is nevertheless bad. The systems labelled 1.3804.164 in Figure [9] and 1.4055.98 in Figure 
[To] show a large and step-like scatter band, the reason of which is, we think, intrinsic variability 
of the secondary componenllfl, as suggested by the fact that the band becomes much narrower at 
primary eclipse but not at secondary eclipse; this interpretation is also suggested by the fact that the 
residuals show an oscillating behavior as a function of phase. In both 1.4055.98 and 36.5943.658 in 
Figure [TT] the scatter band is much larger than the observational error which explains their bad fits. 
These examples show that the samples contain many EBs which could be deserving of more careful 
study which we did not attempt here. The properties of these EBs are summarized in Table [31 of 
the 12 EBs shown 8 have two photometric excursions in both bands and Amply/AmplR < 1.2, 7 are 
found by the decision tree, 5 have s5/s5d < 1.2, and 2 have a counterpart in the OGLE-II sample: 
1.3442.172 (counterpart OGLE050149. 20-691945. 5 in field LMC_SC15) and 1.4055.98 (counterpart 
OGLE050542.06-684732.8 in field LMC_SC13). 

2.6. Ellipsoidal variables in the samples 

Ellipsoidal variability occurs in a close binary system when one (or both) component (s) is 
(are) tidally distorted by the companion. If the binary system is detached, as most systems in our 
samples are, the distorted stars assume the asymmetric, egg like shape of the Roche equipotential 
surface whereas in the case of contact system the shape of the common equipotential surface is more 
reminiscent of a dumbbell. The light curve of an ellipsoidal variable system reveals a continuously 



'The component eclipsed at primary eclipse. 
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Fig. 5.— 



RMS of residuals as a function of median relative magnitude for the LMC sample. 
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Fig. 6. — Examples of LMC EBs light curves, arranged by ascending period; for basic data see 
Table m Left: observed light curves with all data points. The arrows show the baseline as defined 
in Table m Right: observed light curves with outlying points removed, theoretical light curves from 
the fit and residuals. 
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Fig. 7. — Examples of LMC EBs light curves, arranged by ascending period; for basic data see 
Table m Left: observed light curves with all data points. The arrows show the baseline as defined 
in Table m Right: observed light curves with outlying points removed, theoretical light curves from 
the fit and residuals. 
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Fig. 8. — Examples of LMC EBs light curves, arranged by ascending period; for basic data see 
Table m Left: observed light curves with all data points. The arrows show the baseline as defined 
in Table m Right: observed light curves with outlying points removed, theoretical light curves from 
the fit and residuals. 
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Fig. 9. — Examples of LMC EBs light curves, arranged by ascending period; for basic data see 
Table m Left: observed light curves with all data points. The arrows show the baseline as defined 
in Table m Right: observed light curves with outlying points removed, theoretical light curves from 
the fit and residuals. 
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Fig. 10. — Examples of LMC EBs light curves, arranged by ascending period; for basic data see 
Table m Left: observed light curves with all data points. The arrows show the baseline as defined 
in Table m Right: observed light curves with outlying points removed, theoretical light curves from 
the fit and residuals. 
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Fig. 11. — Examples of LMC EBs light curves, arranged by ascending period; for basic data see 
Table m Left: observed light curves with all data points. The arrows show the baseline as defined 
in Table m Right: observed light curves with outlying points removed, theoretical light curves from 
the fit and residuals. 
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Table 4. Basic data for the LMC EBs shown in Figs. [U [TJ and[8l The EEs are arranged by 

ascending period. 



MACHO ID RA(J2000) DEC(J2000) Period^ (d) V baseline'' ''i? baseline'' ''V'-i?'^ Comment 



65.8581.67 


05: 


:33: 


:09, 


,458 


-65: 


:30: 


:29, 


,13 


0.19 


18, 


,93 


18, 


,13 


0.8 


Shortest period in sample 


9.5000.790 


05: 


:11: 


:20, 


,548 


-70: 


:21: 


:00, 


,16 


0.24 


19, 


,70 


19, 


,23 


0.47 


Very short period 


1.3442.172 


05: 


:01: 


:49, 


,005 


-69: 


:19: 


:45, 


,60 


1.02 


17, 


,21 


17, 


,26 


-0.05 


Fairly typical EB 


10.4035.145 


05: 


:05: 


:02, 


,233 


-70: 


:06: 


:13, 


,27 


2.53 


17, 


,14 


17, 


,19 


-0.05 


Fairly typical EB 


68.10843.699 


05: 


:47: 


:10, 


,738 


-67: 


:58: 


:55, 


,74 


3.07 


13, 


,56 


13, 


,99 


-0.43 


Bluest in sample rpj^^ 
Very high eccentric orbit 


62.7240.102 


05: 


:24: 


:55, 


,090 


-66: 


:11: 


:55, 


,28 


4.17 


18, 


,30 


18, 


,30 


0.00 


1.3804.164 


05: 


:03: 


:36, 


,536 


-69: 


:23: 


:32, 


,27 


4.19 


16, 


,88 


16, 


,90 


-0.02 


Algol type 


41.2459.43 


04: 


:55: 


:43, 


,321 


-70: 


:18: 


:00, 


,44 


13.18 


16, 


,74 


16, 


,82 


-0.08 


Highest eccentric orbit 


1.4055.98 


05: 


:05: 


:42, 


,201 


-68: 


:47: 


:33, 


,44 


24.51 


17, 


,32 


17, 


,11 


0.21 


Fairly typical EB 


12.10443.34 


05: 


:44: 


:47, 


,185 


-70: 


:27: 


:26, 


,65 


319.37 


17, 


,35 


16, 


,30 


1.05 


Reddest in sample 


62.6514.2213 


05: 


:20: 


:32, 


,499 


-66: 


:13: 


:17, 


,92 


417.60 


16, 


,88 


16, 


,23 


0.65 


Long Period 


36.5943.658 


05: 


:17: 


:15, 


,478 


-71: 


:57: 


:45, 


,69 


633.70 


16, 


,35 


15, 


,54 


0.81 


Longest Period in sample 



information in Table |4] is also available in its entirety via the link to the machine-readable version above. The EBs 
are arranged by ascending period. Units of right ascension are hours, minutes, and seconds, and units of declination 

are degrees, arcminutes, and arcseconds. 



''Supersmoother provides different periods for the V and R unfolded light curves, but their difference is usually smaller 
than the precision to which we report their values in this table. On line summary tables provide both periods to 5 significant 
digits. 

''The baseline is calculated in the following way. First the outlying points are eliminated by dividing the light curve in 
boxes of 50 data points and eliminating the points which are more than 2 standard deviations away from the mean in each 
box. Then the median of the 10% most luminous points is taken. This value is not the median of the whole light curve that 
is shown in the figures. 

■^Values are quoted to the hundredths of magnitude, typical of MACHO observational uncertainties. 

''This is the difference of the two baselines as defined above, not of the two medians as is the V — R shown in the figures. 
This column is not directly available in the online table but can be deduced by subtracting col. (12) from col. (10). 
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varying profile, with two maxima and two minima per period, with the minima often having different 
depth, whereas the maxima are usually equal. The main reason for this variability is that, as the 
stars rotate, their projected areas on the sky vary, reaching a maximum at the two quadratures 
and a minimum at the two conjunctions; the measured flux thus varies in the sam e way during a 
period. More information on ellipsoidal variables can be found in ([Hilditch l |200l|) A large sample 



of ell ipsoidal variables in the LMC has been released by the OGLE collaboration (ISoszynski et al 



20041'): an ana l ysis o f ellipsoidal variables found in the MACHO database has been published by 
Derekas et al.l ()2006l ). A binary system can present both eclipses and ellipsoidal variability but in 
many cases it may not be possible to clearly recognize an eclipse from visual inspection; this poses 
a problem for the compilation of EB catalogues since the light curves of EBs and non eclipsing 
ellipsoidal variables can be easily confused. 

We looked for possible contamination by ellipsoidal variables in our sample and found ~ 120 
systems exhibiting ellipsoidal variability which we then visually checked more carefully than other 
stars which were clearly E Es. We also a ttempted a less subjec t ive approach by fitting these systems 



with the EBOP program (jEtzel 



1981 



Popper Etzel 



1981 



Nelson fc DavisI Il972l ) following the 



prescriptions of lAlcock et al.l ( 1997al ): however since EBOP is not designed for analyzing such 
distorted systems the final decision about whether or not to include a star exhibiting ellipsoidal 
variability in the sample was taken upon visual inspection. Figure [12] shows two examples of EB 
systems with pronounced ellipsoidal variability, basic data on these systems are given in Table \5\ 



2.7. The Small Magellanic Cloud sample 

The SMC sample comprises 1509 EBs selected via the same techniques as the LMC EBs and 
confirmed by visual inspection; the general considerations of the preceding subsection regarding 
search for variability apply here as well. The sky coverage in the SMC corresponds to MACHO 
fields 206, 207, 208, 211, 212 and 213; field center coordinates for these fields are given in Table 

El 

Magnitudes quoted for the SMC have been obtained by using transformations which dif- 
fer slightly in the zero point from the LMC ones due to the larger exposure times in the SMC 



(jAlcock et al.lll999l ): they are reported in Eq. [2l 



V = yMACHO + 24.97 mag — 0.18(Vmacho ~ -Rmacho) 

R = i?MACHO + 24.73 mag + 0.18(Fmacho - iiMACHo)- (2) 

The SMC search gave 194 duplicates and again we chose the field which had the highest total 
number of observations and the object is identified by that corresponding FTS only. 

Figure [131 shows the histogram of the number of observations in both bands for the EBs in 
the sample. Figure. [Ml shows a logarithmic histogram of the period distribution. Figure [151 
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Fig. 12. — Above: light curves for two long period EBs with strong ellipsoidal variability in the 
LMC sample. Below: light curves of two non eclipsing ellipsoidal variables. 
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Table 5. Basic data for the LMC ellipsoidal variables shown in Figure [T2j The variables are 

arranged by ascending period. 



MACHO ID RA(J2000) DEC(J2000) Period'' (d) F baseline'"' baseline'"' F - i?'^ Eclipsing 



1.3934.140 


05:04:23.977 


-68 


49:21.65 


85.88 


17.68 


17.03 


0.65 


No 


58.6147.42 


05:18:03.677 


-66 


27:20.17 


105.08 


17.78 


17.06 


0.72 


No 


15.10916.25 


05:47:21.831 


-71 


10:46.49 


355.28 


17.04 


16.05 


0.99 


Yes 


14.9588.6 


05:39:28.498 


-71 


00:46.01 


411.04 


14.75 


13.95 


0.80 


Yes 



''Supersmoother provides different periods for the V and R unfolded light curves, but their difference is 
usually smaller than the precision to which we report their values in this table. On line summary tables 
provide both periods to 5 significant digits. 

'^See Table |4] for an explanation of the baseline calculation. 

'^This is the difference of the two baselines, not of the two medians as is the V — R shown in figures. 
''Values are quoted to the hundredths of magnitude, typical of MACHO observational uncertainties. 



Table 6. MACHO field coordinates for the SMC. 



Field ID 


RA(J2000) 


DEC(J2000) 


Date 


206 


1:05:21.70 


-72:26:58.3 


(J2000.0) 


207 


0:57:16.58 


-72:34:57.0 


(J2000.0) 


208 


0:48:03.19 


-72:34:20.9 


(J2000.0) 


211 


0:58:27.40 


-73:04:55.3 


(J2000.0) 


212 


0:49:10.27 


-73:13:32.9 


(J2000.0) 


213 


0:40:18.91 


-73:08:49.5 


(J2000.0) 
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Fig. 13.— 



Histogram of the number of light curve points for both bands for the SMC sample. 
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Fig. 14. — Period histogram for 1509 EBs in the SMC sample. The size of the bins is ~ 1/100 of 
the span of the logarithms of the periods. 
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shows the histograms of the magnitudes for both bands as well as for color. Magnitudes range 
in values from 19 mag to ~ 14 mag both in V and R bands, with a peak around 17 mag. 
The average photometric error for the SMC is again ~ 0.05 mag in both instrumental bands; 
the error as a function of standard magnitude is shown in Figure [T6j Figures [T71 [TE[ [T9\ and [20] 
show some examples of light curves; their properties are summarized in Table [71 The "bump" 
shown by the star labelled 207.16374.39 is probably due to star spots: we were able to roughly 
reproduce it by appropriately choosin g spots on the components and fitting the light curve using 
the PHOEBEH (|Prsa fc Zwitteilbood ) software package. 



3. Color Magnitude Diagram and Color Period Diagram 

3.1. The Large Magellanic Cloud sample 

Figure [21] shows the CMD for the 4634 EBs in the LMC sample; the lower magnitude limit 
is y ~ 21 mag. We estima ted the reddenmg by using the LMC extinction map described in the 
LMC photometric survey o f IZaritskv et al.l (|2004l ). The extinction catalog produced by the survey 
is available for query on hnfl and we retrieved the values of Ay specified in Table [S] based on the 
hot stars only found by the survey (T > 12000K) in a radius of 12' (the maximum allowed) around 
the positions specified in Table [8] which sample the EB position distribution. 

From the values of Table [8] we derive a mean value for Ay of 0.64 mag which we use to 



characterize the average LMC V extinction. We use the reddening vector 



5 from 



Alcock et al.l (jl997bl . and references therein) and find {E{ V — R)) 



E{V-R) 
0.128 mag, more than a factor 



of two and a half larger than the value 0.049 mag found bv I Alcock et al.l (|l997bl ). This is likely due 
to the fact that the re d denin g for the EBs are likely to be along lines of sight toward yo ung, hot stars 
in the IZaritsky et al.l (l2004l ) catalog which are derived to have higher Ay, while the lAlcock et al 



(|l997bl ) Ay was derived from observations of RR Lyrae stars. As the CMD shows the sample is 



made up mostly of bright early type stars: from the rang e in magnitudes and assuming an LMC 



distance modulus of 18.5 mag (Ivan der Marel et al 



component is of spectral type B or A (lAlcock et al 



2002) we see that in most cases at least one 



2000bl : ICoxll2000l ). We employ the term young 



star region to describe the main feature on the left part of the CMD, rather than Main Sequence 
because our sample contains some bright and short lived stars that may not be burning hydrogen 
in their core while still being on the blue part of the CMD; likewise we employ the term evolved 
star region to indicate the feature on the red part of the CMD. We define the young star region as 
V — R< 0.2 mag and the evolved star region as V — R > 0.2 mag: with these definitions we find 
3760 EBs in the young star region and 874 EBs in the evolved star region. We used a simple cut to 
separate the young star region from the evolved star region rather than a more precise one because 



^" ^http: //phoebe .f iz .uni-lj . si/| 

^ ^http: //ngala. as . arizona. edu/dennis/lmcext .html| 
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Fig. 15. — V and R magnitudes and V — R histograms for 1508 EBs in the SMC sample. The bin 
size is 0.1 mag for the V and R histograms and 0.01 mag for the V—R one. 
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Fig. 16.— 



Errors as a function of median magnitude for the SMC sample. 
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Fig. 17. — Examples of light curves of EBs in the SMC sample; for basic data see Table [71 The 
arrows show the baseline as defined in Table [71 The EBs are arranged by ascending period. 
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Fig. 18. — Examples of light curves of EBs in the SMC sample; for basic data see Table [71 The 
arrows show the baseline as defined in Table [71 The EBs are arranged by ascending period. 
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Fig. 19. — Examples of light curves of EBs in the SMC sample; for basic data see Table [71 The 
arrows show the baseline as defined in Table [71 The EBs are arranged by ascending period. 
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Fig. 20. — Examples of light curves of EBs in the SMC sample; for basic data see Table [71 The 
arrows show the baseline as defined in Table [71 The EBs are arranged by ascending period. 
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Table 7. Basic data for the SMC EBs shown in Figs. [HI [HI 19] and [20l The EEs are arranged 

by ascending period. 



MACHO ID 


RA(J2000) 


DEC(J2000) 


Period^ (d) 


V baseline'' 


''R baseline'' 




Comment 


211.16529.5 


00:59:31.368 


-73:26:56.04 


0.28 


16.00 


15.43 


0.57 


Shortest period in sample 


207.16652.1084 


01:00:43.656 


-72:51:51.48 


0.36 


19.67 


19.59 


0.08 


Very short period 


206.16883.214 


01:04:25.272 


-72:37:48.36 


1.14 


18.16 


18.18 


-0.02 


Fairly typical EB 


207.16656.93 


01:00:56.345 


-72:36:44.40 


1.23 


18.00 


17.96 


0.04 


Highly eccentric orbit 


206.16717.285 


01:01:41.400 


-72:19:25.68 


1.65 


18.35 


18.32 


0.03 


Fairly typical EB 


208.15740.63 


00:46:13.949 


-72:52:37.03 


1.74 


16.13 


16.24 


-0.11 


Bluest in sample 


212.15903.2269 


00:49:18.192 


-73:21:55.44 


2.42 


17.51 


17.51 


0.00 


Highly eccentric ^iJsit 


207.16315.289 


00:55:48.168 


-72:29:33.36 


3.34 


17.93 


17.94 


-0.01 


Highly eccentric orbit 


211.16195.61 


00:53:59.729 


-72:56:56.13 


4.73 


16.44 


16.43 


0.01 


Fairly typical EB 


208.15912.323 


00:49:10.440 


-72:46:37.56 


120.51 


18.84 


18.35 


0.49 


Highly eccentric orbit 


208.16.58 


00:49:28.392 


-72:49:40.80 


137.89 


17.17 


16.50 


0.67 


Reddest in sample 


207.16374.39'* 


00:56:25.872 


-72:22:15.96 


186.34 


16.30 


16.30 


0.00 


Long Period 


212.15673.13 


00:45:46.824 


-73:31:32.52 


200.27 


16.00 


15.35 


0.65 


Long period 


211.16418.53 


00:57:5.304 


-73:15:10.44 


234.64 


17.23 


16.82 


0.41 


Highly eccentric orbit 


206.17005.6 


01:06:10.224 


-72:06:24.48 


371.89 


15.69 


15.16 


0.53 


Highly eccentric orbit 


211.16310.174 


00:55:30.024 


-72:52:49.44 


1559.81 


18.23 


17.61 


0.62 


Longest Period in sample 



information in Table [7] is also available in its entirety via the link to the machine-readable version above. The EBs 
are arranged by ascending period. Units of right ascension are hours, minutes, and seconds, and units of declination 

are degrees, arcminutes, and arcseconds. 



''Supersmoother provides different periods for V and R unfolded light curves, but their difference is usually smaller than the 
precision to which we report their values in this table. On line summary tables provide both periods to 5 significant digits. 

''The baseline is calculated in the following way. First outlying points have been eliminated by dividing the light curve in 
boxes of ~ 50 data points and eliminating the points which were more than 2 standard deviations away from the mean in each 
box. Then the median of the 10% most luminous points was taken. This value is not the median of the whole light curve that 
is shown in the figures. 

'^Values are quoted to the hundredths of magnitude, typical of MACHO observational uncertainties. 

''This is the difference of the two baselines as defined above, not of the two medians as is the V — R shown in the figures. This 
column is not directly available in the online table but can be deduced by subtracting col. (12) from col. (10) 

■^There is a curious "bump" in the light curve of this long period EB (Figure I19p that suggests further investigation. 
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our cut can be easily seen both in the CMD and in the Color Period Diagrams. An interesting 
feature of the CMD is the lack of a clear gap between the young star region and the evolved star 
region; there is instead a fairly continuous transition, with a higher number of systems filling the 
Hertzsprung Gap that would be expected from CMDs of single stars. This may indicate that these 
systems are composed of a more massive and hence more evolved and redder star and a less massive, 
less evolved, bluer one. 



The Color Period Diagram is shown in the left panel of Figure [22l this diagram clearly shows the 
young EBs with periods P < 20d and a second population (333 objects) of long period, evolved EBs 
with periods P > 20d and V — R > 0.2 mag. Several interesting features emerge in Figure [22t there 
is a paucity of long period objects on the young star region of the CMD, and a corresponding lack 
of short period objects with very red colors. Furthermore, the red {V — R> 0.2 mag) population 
shows a positive correlation between period and color. There are virtually no long period, blue 
objects or short period, red objects in this group. In contrast, the young star region shows no 
such correlation. This structure in the Color Period Diagram is a consequence of (i) Kepler's 
Third Law, (ii) the probability that an EB is favorably oriented in space to allow eclipses to be 
detected (Prob = {Ri + R2)/a) where Ri and R2 are the radii of the primary and secondary stars, 
respectively and a is the semi major axis, and (iii) that objects usually evolve in this diagram at 
constant period, from blue to red. Long period binary stars have large semi-major axes. When 
both stars are on the young star region, their relatively small radii yield a relatively low probability 
that they will eclipse when seen from our vantage point. When one of the pair evolves away from 
the young star region, one of these radii (typically Ri) will increase. The consequence of this is 
an increase in the probability that eclipses will be detected. This accounts for the presence of red, 
long period stars and the absence of corresponding young progenitors. The situation is different for 
short period systems. These are relatively likely to be detected because of their small semi-major 
axes, and are prominent on the blue side. As one of these stars evolves and expands rapidly, it may 
engulf the companion and enter a stage o f common envelo pe evolution in which the expanding star 
overflows the second Lagrangian point L2 (|Paczvnskilll976); leading to the disappea r ance of eclipses . 



Common envelope system differ fro m contact binaries (jKallrath &: Milond Il999l : iHilditch 1 12001 



Shore. Livio &: van den Heuvell Il994lri in which two young stars overflow their first Lagrangian 
point (Li) and their Roche equipotential surface assumes a dumbbell shape. Contact systems 
usually show ellipsoidal variation and, if t he orbital inclination is large enough, also eclipses; EBs 
of the W UMa type belong to this category. IShore. Livio &: van den Heuvell (jl994l ) and llben &: Livio 
(|l993l ) provide more information on common envelope binaries. The correlation between period and 
color among red objects reflects the general correlation between radius and color for the evolved 
partner. 



i^Called "over contact" in l|Kallrath fc Milond [l999l ) 
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0.5 1 

V-R 

Fig. 21. — CMD for 4634 EBs in the LMC sample. The sample is made up mostly by young 
luminous stars (with at least one component of spectral type B or A). A fairly high number of EBs 
(~ 19%) are evolved; there are no real "gaps" between the young star region and the evolved star 
region. The vertical line shows the cut used to separate the young star region from the evolved 
star one. The reddening vector is ^^y^ = 5; the adopted values of {Ay) and {E{V — R)) are 
0.64 mag and 0.128 mag. 
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Table 8. Extinction values in the LMC. 



HA(J2000) 


Di!jC(J20Q0) 


Av 


06:06:00 


-69:05:00 


0.48 


06:06:00 


-72:43:00 


0.97 


05:40:00 


-65:30:00 


0.75 


05:40:00 


-69:05:00 


0.80 


05:40:00 


-72:30:00 


0.78 


05:20:00 


-65:30:00 


0.53 


05:20:00 


-69:05:00 


0.50 


05:20:00 


-72:30:00 


0.56 


05:00:00 


-65:30:00 


0.50 


05:00:00 


-69:05:00 


0.46 


05:00:00 


-72:43:00 


0.45 


04:40:00 


-69:05:00 


0.62 


04:40:00 


-72:30:00 


0.97 




Fig. 22. — Left Panel: Color Period Diagram for 4634 EBs in the LMC. The gray filled triangles 
represent the foreground population. Right Panel: Color Period Diagram for 1508 EBs in the SMC. 
The figure shows the higher fraction of long period EBs that belong to the young star region in the 
SMC than in the LMC. 
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3.2. Foreground objects 

The left panel of Figure [22] reveals a population of ~ 63 EBs with low periods {P < 2d) and 
high color {V — R> 0.3 mag). These objects are probably foreground galactic EBs composed of 
late type stars. This interpretation is suggested by several factors. First, due to the large angular 



extent of the LMC, there is foreground contamination in the LMC MACHO fields (lAlcock et al 



2000b!); in particular the feature marked "H" in the CMD of their Figure 1 indicates foreground 
galactic disk stars and is centered at V — R^ 0.5 mag as is our presumptive foreground population. 
Second, both the short period of these EBs, and the shape of their light curves which are either 
detached or mildly distorted, strongly suggests that the stars making up this population are small, 
late type stars; this is further borne out by their color, again typical of a solar like star; Finally, 
the CMD of this population, shown in Figure [23l clearly shows what appears to be a turnoff feature 
at V — R^ 0.4 mag, with few evolved objects {V — R> 0.7 mag). It is interesting to note that the 
overall shape of this population in the Color Period Diagram shows, on a smaller scale, the same 
features of the LMC diagram; a Main SequencJ^ is clearly visible in Figure [23] and the evolved 
EBs show the same Color Period correlation of their LMC counterparts in Figure [22] 



3.3. The Small Magellanic Cloud sample 

Figure [M] shows the CMD for 1508 EBs in the SMC sample out of the 1509 in the sample; 
one EB which has valid data only in the iiMACHO band is not shown because it was not possible 
to determine the standard magnitudes via Eq. [2] The general remarks made for the LMC CMD 
apply here as well and we used a reddening vector with the same inclination as the LMC. The figure 
clearly shows the young star region which is composed o f 1412 EBs (94%) wh ereas there are just 



96 evolved EBs. We estimate the SMC reddening from IZaritskv et al.l (|2002l ): from their Figure 
19 we infer a mean Ay ~ 0.3 mag for their hotter SMC population, relevant to our sample which 
is composed mostly of early type hot stars on the young star region. We use the same reddening 
vector as the LMC, ^^^y^ji-j ~ ^ mean E{V — R) = 0.06 mag. The Color Period 

Diagram is shown in the right panel of Figure [22l the general considerations made for the LMC 
Color Period Diagram apply here as well. The diagram shows 5 EBs with low period (P < Id) and 
red color (V — R^ 0.6 mag) which are probably foreground objects. The much smaller number of 
foreground objects in the SMC is probably due both to its smaller angular size and to its higher 
galactic latitude {I ~ —44° compared to I ~ —33° for the LMC). We can further test the hypothesis 
that these two short period, red populations in the LMC and SMC are due to foreground objects 
by comparing the ratio of their numbers to the ratio of the areas of the LMC and SMC, which we 
can estimate from Figures [28l and [37] if these objects are indeed foreground these two ratios should 
be roughly equal. From the figures we can estimate the sky area of the LMC as ~ 110n° and the 



We employ this term since these foreground stars are probably not very massive and therefore are in their core 
hydrogen burning phase. 
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Fig. 23. — CMD of the LMC sample (left panel) and of the foreground population (right panel,filled 
hexagons). 
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Fig. 24.— CMD for 1508 EBs in the SMC sample. The reddening vector is ^^y^Lji^ = 5; the 
adopted values of {Ay) and {E{V — R)) are 0.3 mag and 0.06 mag. 
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sky area of the SMC as ~ 10n°; their ratio is thus similar to the ratio of the numbers of EBs in 
the two populations as expected. 



4. Comparison between the Large Magellanic Cloud and the Small Magellanic 

Cloud samples 

Although the basic features of the CMD and the Color Period Diagram are the same for LMC 
and SMC, there are some differences, shown by Figure \n\ and Table 

The fraction of blue {V-R < 0.2 mag) EBs is higher in the SMC than in the LMC. More 
striking, the fraction of blue long period {V — R< —0.2 mag, P > 20d) EBs is much higher in the 
SMC than in the LMC. To further investigate the differences between the two samples we carried 
out a Kolmogorov-Smirnov (KS) test on the distributions of the absolute magnitudes My and Mr, 
their difference My — Mr , and the periods P of the two samples. We adopted a distance inodulu s 



of of 18.88 for the SMC (|Dolphin et al.ll2n0lh and 18.5 for the LMC (|van der Marel et al.l 120021 ). 

Magnitudes and colors were dereddened using (Ay) = 0.64 mag, {E{V — R)) = 0.128 mag for the 
LMC and {Ay) = 0.3 mag, {E{V — R)) = 0.06 mag for the SMC before subtracting the distance 
moduli. The Empirical Cumulative Distribution Functions (ECDFs) of these quantities are shown 
in Figure [25l A KS test confirms that the distributions of P and My — Mpi are different at > 99.9% 
confidence level; for the distributions of My and Mji the KS test gives a probability of ~ 99.4% 
and ~ 98.5% respectively for them being different. The My — Mr plot shows that EBs in the 
SMC tend to be bluer, not surprising given that a higher percentage of objects belong to the young 
star CMD region in the SMC than in the LMC. The period distributions show that the SMC EBs 
have on average shorter periods, again not surprising given the much higher percentage of evolved 
systems in the LMC than in the SMC and the fact that evolved systems have on average higher 
periods than the young systems (as shown by the Color Period Diagrams in Figure [22]) . 

Table 9. Summary of long period EBs in the SMC and LMC. 



Galaxy 


Total 


Young stars ^ 


Evolved ^ 


Long Period 


Long Period young stars 


^ Long Period evolved stars ^ 


LMC 


4634 


3760(81%) 


874(19%) 


356 


23(6%) 


333(94%) 


SMC 


1508 


1412(94%)'^ 


96(6%)'' 


75 


23(31%) 


52(69%) 



^Defined as V-R < 0.2 mag. 
''Defined as V'-i? > 0.2 mag. 
'^Defined as P > 20d. 

''One SMC EB has no valid Vmacho data, hence the sum of the young star and evolved star numbers for the SMC is 1508. 




Fig. 25.— ECDFs for the distribution of My (above left), Mr (above right), My - Mr (Below 
left) and P (below right) for the SMC (continuous line) and the LMC (dashed line). 
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5. Cross correlation with the OGLE-II samples 

The EBs in our samples have been cross correlated with the corresponding OGLE-II samples. 
Stars in the samples were identified if their right ascension (RA) and declination (DEC) differed 
by less than 27.2" and if their periods differed by < 1%. We used a very large search radius to be 
conservative. The astrometric precision for both surveys is typically ~ 1", but a few stars which 
had much larger differences in RA and/or DEC turned out to be matches upon inspection of their 
periods: in particular we found in the LMC 32 matches with a position difference bigger than 10" 
and 2 matches with a position difference bigger than 20". However most of the matches were within 
narrower radii: for the LMC roughly half of the matches (534 out of 1236) were found within of 
~ 2", compatible with the astrometric precision of both MACHO and OGLE surveys; almost all 
of them (1019) were within ~ 4". We tested the robustness of our method of finding matches by 
investigating the probability for two periodic objects with a period difference of < 1% to be within 
27.2" of each other. To do this we selected a random sample of 5000 objects out of the ~ 66000 
periodic ones found by MACHO in the LMC and we counted the frequency of pairs of objects 
with both periods from the red and the blue lightcurves (as found by Supersmoother) differing 
by < 1% and positions within 27.2". Most of the matches we found were due to the same object 
being observed in different tiles and only in one case did we find a possibly genuine match; we thus 
conclude that the probability of two objects being erroneously classified as a match is ~ 0.02% and 
therefore our method of finding matches is robust. 

5.1. The Large Magellanic Cloud sample 

Our search produced 1236 matches in the LMC. The MACHO and OGLE-II periods agree to 
high accuracy as shown in Figure [26l usually much better than the 1% cut we imposed. Histograms 
of position differences are shown by Figure [271 Both panels show the entire span of the differences; 
the differences in Right Ascension range from —20" to -|-20" but the left panel shows this 
range multiplied by the cosine of the declination (~ —70°) which gives a range from ~ —7" to 
~ -|-6". The sky coverage of the two surveys was different: OGLE-II, from which the sample was 
derived, covered about 4.5n° in the central region of the LMC whereas the sky coverage of MACHO 
was larger. Figure [28] shows the positions of the EBs in both catalogues and Figure [29] shows the 
corresponding MACHO field number; the fields at the center of the LMC are drawn with continuous 
lines, the ones at the periphery with dashed lines. 

5.2. Comparison between the center and the periphery of the LMC 

In view of the different sky coverage of the MACHO and OGLE-II surveys, it is interesting 
to analyze separately the stars in MACHO fields covering the center of the LMC (which roughly 
correspond to the OGLE-II sky coverage) and the stars in the MACHO fields at the periphery. 
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Fig. 26.— Percentage difference for MACHO vs. OGLE-II period for the 1236 OGLE-II matches 
in the LMC sample. 
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Figure. [29] shows that the fields in the center are 1, 2, 3, 5, 6, 7, 8, 10, 11, 12, 13, 14, 15, 18, 
19, 47, 77, 78, 79, 80, 81 and 82. There are 2620 EBs in the center and 2014 in the periphery. 
Figure [30] shows the CMD and the Color Period Diagram for the center and the periphery of the 
LMC respectively; Figure [31] shows the histograms for the magnitudes, the color and the period. 
The figures reveal several differences between the two samples; to check these we performed a 
KS test for V, V—R and P and found that at high confidence level (> 99.9%) the distributions 
are different. These differences are statistically significant but not large enough to be considered 
astrophysically important. In particular the difference in period distribution could be attributed 
to the lower sampling in the outer fields, making it less likely to detect the longer period EBs. The 
sampling might also be the cause of the periphery having a higher percentage of bright EBs, since 
these are easier to find with fewer epochs. 



5.3. Discussion of OGLE-II MACHO comparison 

We finally investigated why we did not find more matches with OGLE-II. The most important 
reason, we think, is the fact that the techniques employed i n assembling the samples are di fferent: 



the OGLE team built their samples via neural networks (IWyrzykowski et al.l l2003l . [200J). The 
two surveys have roughly comparable limiting magnitudes, V ^ 21.5 mag; nevertheless we checked 
how the performance of the two surveys varied with magnitude. We compared the distribution of 
the V magnitudes for the 2620 EBs in the central region of the LMC in our sample with those of 
1198 matches and the 1327 EBs from OGLE-II without MACHO counterparts (these two numbers 
do not add up to 2580, the size of the OGLE-II sample, because for some EBs the V magnitude 
was not reported). The histograms of these three distributions are shown in Figure [32] The 
figure shows that OGLE-II EBs with MACHO counterparts, peaking at ^ ~ 18 mag like the 
MACHO sample, are on average brighter than the ones without MACHO counterparts, which peak 
at y ~ 19.5 mag. The shape of the V distribution of OGLE-II EBs with MACHO counterparts 
much more closely resembles the MACHO V distribution; a KS test gives a probability of the two 
distributions being the same of ~ 46%. The distributions of MACHO V and OGLE-II V without 
MACHO counterparts, as well as those of these two OGLE-II populations are, on the other hand, 
shown to be different at > 99.9% confidence level. Both distributions vanish at y ~ 20.5 mag, 
showing that their limiting magnitudes are comparable. 

We then studied the distributions of the periods: Figure [33] shows that the OGLE-II EBs 
without MACHO counterparts have on average longer periods than MACHO EBs in the center of 
the LMC and than OGLE-II EBs with a MACHO counterpart and the difference is statistically 
significant in both cases; the period distributions of the MACHO EB of the LMC center and of the 
MACHO OGLE-II matches are statistically different as well. Therefore, OGLE-II finds a higher 
proportion of fainter objects than MACHO and this does play a role in not finding an higher number 
of OGLE-II counterparts to our sample. We finally studied the distribution of the MACHO-OGLE 
matches as a function of V. We first counted the number of OGLE-II LMC EBs in the MACHO 
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fields, finding 2517 of tliem, out of a total of 2580; of these 2517 EBs, 1236 were the matches 
described above and 1281 did not have a MACHO counterpart (this last number is smaller than 
1327, the total number of OGLE EBs without MACHO counterpart, because we are now only 
considering OGLE EBs in MACHO fields). We then studied the distribution, as a function of V, of 
the OGLE EBs in the MACHO fields that both had and did not have a MACHO counterpart and 
for which the V magnitude was reported: there were 1198 of the former and 1267 of the latter. We 
finally performed the inverse calculation, by first counting the number of MACHO LMC EBs in the 
OGLE-H fields and finding 1551 of them; of these 1225 had an OGLE-II counterpart and 326 did 
not; we studied the distribution of both these populations as a function of V. The OGLE-II field 
boundaries were estimated by taking the coordinates of the most extreme EBs in each OGLE-II 
field. These findings are summarized in Table [TOj 

The distributions of both the OGLE-II EBs with MACHO counterpart and of the MACHO 
EBs with OGLE-II counterparts are shown in Figure [34l the figure shows the fraction of matches 
in magnitude bins of 1 mag; the bin centers range from V = 20.5 mag to V = 12.5 mag. The 
error bars are estimated by assuming that the matches in each magnitude bin follow a binomial 
distribution with probability p = x/Nj, where x is the number of matches in each magnitude bin 
and Njj is the total number of EBs; the error in the expected fraction of matches is then given by 
EqO in both distributions of Figure. [3l]the error bar for the brightest magnitude bin is not shown 
since in both cases there is only one match, rendering Eq. [3] meaningless. The figure shows that 
the fraction of matches increases for brighter magnitudes as expected; the fall at V = 13.5 mag in 
the distribution of OGLE matches is probably due to small number statistic as evidenced by the 
large error bar. 




(3) 



Table 10. MACHO OGLE-II matches. 



OGLE-II LMC EBs 2580 

OGLE-II LMC EBs in MACHO fields 2517 

OGLE-MACHO matches 1236 

OGLE-II LMC EBs without MACHO counterpart 1281 

OGLE- MACHO matches with reported V 1198 
OGLE-II LMC EBs without MACHO counterpart with reported V 1267 

MACHO EBs in OGLE fields 1551 

MACHO EBs in OGLE fields with OGLE counterpart 1225 

MACHO EBs in OGLE fields without OGLE counterpart 326 
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5.4. The Small Magellanic Cloud sample 

The same general considerations apply to the SMC sample: the search, performed with the 
same criteria as the LMC, produced 698 matches. Figure [35] shows the percentage difference of the 
MACHO and OGLE-II periods vs MACHO period for the matches; Figure [36] shows the histogram 
of the differences in RA and DEC. Unlike the LMC, the sky coverage of the two surveys was 
approximately the same. Figure [37] shows the positions of the MACHO and OGLE EEs on the 
sky. We again investigated why we did not find more matches with OGLE-H. Looking at Figure 
[37] it is evident that one of the reasons is the somewhat different sky coverage of the two surveys 
(MACHO Fields 207, 208 and 211 are only partially covered by OGLE), but we also looked for other 
possible explanations. Again, the most likely explanation is the ways in which the samples were 
assembled, but we also considered the differences in the distributions of magnitudes and periods. 
As for the LMC, we compared the distribution of the V magnitudes for the 1508 EEs in our 
sample which have a valid V with those of 650 matches and the 666 EEs from OGLE-II without 
MACHO counterparts (again these two numbers do not add up to 1351, the size of the OGLE- 
II sample, because for some EEs the V magnitude was not reported). The histograms of these 
three distributions are shown in Figure [38l The figure shows that OGLE-II EEs with MACHO 
counterparts, peaking at y ~ 17.5 mag like the MACHO sample, are on average brighter than 
the ones without MACHO counterparts, which is more spread out and roughly constant between 
18.5 mag < V < 16.5 mag. A KS test shows that these three distributions are different at > 99% 
confidence level. The behavior at high magnitudes is different for MACHO and OGLE-II: the 
figure suggests a magnitude limit of y ~ 21 mag for MACHO and y ~ 20 mag for OGLE-II; for 

V <19 mag MACHO finds many more EEs than OGLE-II. 

We then studied the distributions of the periods: Figure [39] shows that OGLE-II EEs with 
MACHO counterparts have periods that cluster more in the Id < P < lOd range, as do the MACHO 
EEs and a KS test gives a probability ~ 43% for the two distribution to be the same. On the other 
hand OGLE-II EEs without MACHO counterparts have a larger spread of period values, with more 
objects having P > lOd and with small "bumps" in the distributions at ~ 20d and lOOd and a 
KS test shows it to be different from both MACHO and OGLE-II MACHO matches at > 99.9% 
confidence level. The V distributions of both the OGLE-II EEs without MACHO counterpart 
and of the MACHO EEs without OGLE-II counterparts are shown if Figure [^ the figure shows 
the fraction of matches in magnitude bins of 1 mag with centers ranging from V = 20.5 mag to 

V = 14.5 mag; the error bars are estimated again by Eq. [3] 

We conclude that differences in sky coverage and in techniques used in assembling the samples, 
as well as different magnitude limits and in general different behavior at high magnitude of the two 
surveys all play a role in not finding an higher number of OGLE-II counterparts to our sample; 
the fact that the SMC was less observed than the LMC by MACHO may also explain why we find 
fewer long period objects, wh ile the fact that the exposure time for the SMC was double that of 



the LMC (lAlcock et al.l[l999l ) may explain why we find many faint objects. 
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6. The data on line 

The data presented in this paper can be accessed on hne at the Astronomical Journal websitJ^ 
and are mirrored at the Harvard University Initiative in Innovative Computing (IIC) /Time Series 
CenterH At both sites the data consist of a summary table for each Cloud and light curves for all 
the EBs in the samples. Light curve files contain the unfolded data in MACHO magnitudes; these 
same files can also be retrieved from the MACHO website. Finally for the LMC, the input and 
output files used in the JKTEBOP fits are provided. 
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Fig. 27. — Left Panel: histogram of the differences between Right Ascensions for 1236 MACHO 
OGLE-II matches in the LMC. Right Paneh histogram of the differences between dechnations. The 
bin size is equal to 1/30 of the range of the differences in both cases. 
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Fig. 28. — Cross correlation between MACHO and OGLE-II LMC samples: points represent MA- 
CHO stars, crosses OGLE stars and gray empty boxes the matches. 
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Fig. 29.— 



MACHO LMC field numbers. 
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Fig. 30.— Upper Left Panel: CMD for 2620 EBs in the center of the LMC. Upper Right paneh 
CMD for 2014 EBs at the periphery of the LMC. The CMDs suggest a more continuous transition 
from the young star region to the evolved star region in the center than in the periphery, especially 
for F < 19 mag. 

Lower panels: Color Period Diagrams for the same populations. Lower Left Panel: center. Lower 
Right panel: periphery. The Color Period Diagrams reveal the presence of a long period (20 — lOOd), 
relatively unevolved {V — R'^ 0.2 mag) population in the center but not in the periphery. 
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Fig. 31. — Upper Left Panel: Period histogram for 2620 EBs in the center of the LMC. Upper 
Right panel: period histogram for 2014 EBs at the periphery of the LMC. The size of the bins is 
~ 1/100 of the span of the logarithms of the periods. 

Lower panels: Magnitude and Color histograms. Lower Left Panel: center. Lower Right panel: 
periphery. The bin size is 0.1 mag for the V and R histograms and 0.01 mag for the V — R one. 
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Fig. 32. — V distribution for MACHO 2620 EBs in the central region of the LMC (continuous 
hne), 1198 OGLE-II MACHO matches (dashed hne), and 1327 OGLE-II EBs without MACHO 
counterparts (long dashed line). The figure shows that OGLE-II EBs with MACHO counterparts 
are on average brighter than the ones without, and the shape of their V distribution more closely 
resembles the MACHO V distribution. This is confirmed by a KS test. 
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Fig. 33. — Period distribution for 2620 MACHO EBs in the central region of the LMC (continu- 
ous Hne), 1198 OGLE-MACHO matches (dashed hne), and 1327 OGLE-H EBs without MACHO 
counterparts (long dashed line). The figure shows that OGLE-H EBs with MACHO counterparts 
have on average shorter periods than the ones without. 
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Fig. 34. — Continuous line: V distribution of the fraction of matches for the OGLE-II LMC EBs 
in MACHO fields with expected error bars. The magnitude bins are 1 mag wide and their centers 
range from V = 20.5 mag to V = 12.5 mag. Dashed line: V distribution of the fraction of matches 
for the MACHO LMC EBs in OGLE-II fields. 
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Fig. 35. — Percentage difference for MACHO vs. OGLE-II period for the 698 OGLE-II matches in 
the SMC sample. 
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Fig. 36. — Left Panel: histogram of the differences between Right Ascensions and for 698 OGLE 
matches in the SMC sample. Right Panel: histogram of the differences between declinations. 
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Fig. 37. — Cross correlation between MACHO and OGLE-II SMC samples: points represent MA- 
CHO stars, crosses OGLE stars and gray empty boxes the matches. 
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Fig. 38.— Continuous line: V distribution for 1508 MACHO EBs in the SMC. Dashed line: V dis- 
tribution for 650 OGLE-MACHO matches. Long dashed line: V distribution for 666 OGLE-H EBs 
without MACHO counterparts. The figure shows that OGLE-H EBs with MACHO counterparts 
are on average brighter than the ones without, and the shape of their V distribution more closely 
resembles the MACHO V distribution. 
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Fig. 39.— Continuous line: Period distribution for 1508 MACHO EBs in the SMC. Dashed hne: 
Period distribution for 650 OGLE-MACHO matches. Long dashed Une: Period distribution for 666 
OGLE-H EBs without MACHO counterparts. The figure shows that both MACHO and OGLE-H 
EBs with MACHO counterparts have periods that cluster more in the Id < P < lOd range. OGLE 
EBs without MACHO counterparts have a larger spread in period and smaller "bumps" in the 
distribution at ~ 20d and lOOd. 
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Fig. 40. — Continuous line: V distribution of the fraction of matches of matches for the OGLE-II 
SMC EBs in MACHO fields with expected error bars. The magnitude bins are 1 mag wide and 
their centers range from V = 20.5 mag to F = 13.5 mag. Dashed line: V distribution of the 
fraction of matches for the MACHO SMC EBs in OGLE-H fields. 



